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Abstract. We generalize the results of [Comm. Math. Phys. 299 (2010), 825-866] (hidden 
Grassmann structure IV) to the case of excited states of the transfer matrix of the six- 
vertex model acting in the so-called Matsubara direction. We establish an equivalence 
between a scaling limit of the partition function of the six- vertex model on a cylinder with 
quasi-local operators inserted and special boundary conditions, corresponding to particle- 
hole excitations, on the one hand, and certain three-point correlation functions of conformal 
field theory (CFT) on the other hand. As in hidden Grassmann structure IV, the fermionic 
basis developed in previous papers and its conformal limit are used for a description of the 
quasi-local operators. In paper IV we claimed that in the conformal limit the fermionic 
creation operators generate a basis equivalent to the basis of the descendant states in the 
conformal field theory modulo integrals of motion suggested by A. Zamolodchikov (1987). 
Here we argue that, in order to completely determine the transformation between the above 
fermionic basis and the basis of descendants in the CFT, we need to involve excitations. On 
the side of the lattice model we use the excited-state TBA approach. We consider in detail 
the case of the descendant at level 8. 

Key words: integrable models; six vertex model; XXZ spin chain; fermionic basis, thermo- 
dynamic Bethe ansatz; excited states; conformal field theory; Virasoro algebra 
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1 Introduction 

Much progress was made in the understanding of the connection between the one-dimensional 
XXZ spin chain and two-dimensional quantum field theories (QFT). Many different aspects 
of this connection were studied in the literature. In the present paper we will touch only 
some particular aspect related to the hidden fermionic structure of the XXZ model and the 
corresponding continuum model - conformal field theory (CFT), [IJ [2j [3l U] . 

The continuum model can be studied through the scaling limit of the six-vertex model com- 
pactified on an infinite cylinder of radius R. The corresponding direction around the cylinder 
is sometimes called Matsubara direction. It is well known that the scaling limit of the homoge- 
neous critical six- vertex model is related to CFT, while introducing inhomogeneities in a special 
way leads to the sine-Gordon model (sG) which has a mass gap, [5]. An important step forward 
was done by Bazhanov, Lukyanov and Zamolodchikov in papers [U [7J [8] . They obtained an 
integrable structure of CFT by constructing a monodromy matrix with the quantum space re- 
lated to the chiral bosonic field and the Heisenberg algebra. On the other hand this monodromy 

*This paper is a contribution to the Special Issue "Relationship of Orthogonal Polynomials and Spe- 
cial Functions with Quantum Groups and Integrable Systems". The full collection is available at 
|http://www.emis.de/journals/SIGMA/OPSF.html| 
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matrix satisfies the Yang-Baxter equation with the i?-matrix of the six- vertex model. The corre- 
sponding transfer matrix fulfills Baxter's TQ-relation [U] and generates the integrals of motion. 
These integrals of motion are special combinations of the Virasoro generators. Originally they 
were introduced by Alexander Zamolodchikov in [10] . Zamolodchikov observed that the inte- 
grals of motion generate that part of the Virasoro algebra that survives under the integrable 
^1,3-perturbation of the CFT. 

In order to state the full equivalence of the six-vertex model in the scaling limit and the 
CFT one needs to compare all possible correlation functions. This problem is far from being 
solved completely. We believe that it can be helpful to use a hidden fermionic structure of the 
spin-| XXZ chain, [H [2]. The key idea is to consider a fermionic basis generated by means 
of certain creation operators t*, b*, c*. These operators act on a space of quasi-local oper- 
ators W( a ). Any operator in this space can be represented in the form q 2aS ( )Q with some 
operator which acts on a chain segment of arbitrary finite length. The operator 

q 2aS(0) wag 

called "primary field" because it fulfills some properties similar to the properties of certain 
primary field in CFT. The fermionic basis is constructed inductively. Its completeness was 
shown in [IT]. An important theorem proved by Jimbo, Miwa and Smirnov in the paper [3] 
allows one to reduce any correlation function expressed through the fermionic basis to determi- 
nants. An interesting feature of the above construction is that it is algebraic in the sense that 
the fermionic operators are constructed by means of the representation theory of the quantum 
group U q {s\2)- They are independent of any physical data like the magnetic field, the temper- 
ature, the boundary conditions etc. For example, the temperature can be incorporated via the 
Suzuki-Trotter formalism |12] by taking inhomogeneity parameters in the Matsubara direction 
in a special way and then performing the so-called Trotter limit, when the number of sites in 
the Matsubara direction n — > 00. One can even keep the number n finite and consider the 
case of arbitrary inhomogeneity parameters in this direction. The fermionic basis will not be 
affected [3]. There are only two transcendental functions p and co which "absorb" the whole 
physical information and appear in the determinants in analogy with the situation with free 
fermions. In [3] these functions were represented in terms of the differential of the second kind 
in the theory of deformed Abelian integrals for finite n. These functions can be also obtained 
within the TBA approach, [13] which makes it possible to take the Trotter limit n — > 00. In 
the recent papers [14} [T5] the fermionic basis was used in order to study the one-point func- 
tions of the sine-Gordon model on a cylinder. In this case the transcendental functions must be 
modified. 

Coming back to the similarity with the CFT, let us emphasize that the construction of the 
fermionic basis bears some similarity with the construction of the descendants in CFT via the 
action of the Virasoro generators on the primary fields. In the paper [1] we tried to make this 
similarity more explicit. Namely, we related the main building blocks of CFT - the three-point 
correlators to the scaling limit of a special partition functions of the six- vertex model constructed 
with the help of the fermionic operators. More concretely, we believe that the following con- 
jecture is true. There exist scaling limits t* , (3*, 7* of the operators t*, b*, c* and scaling 
limits p sc and w sc of the functions p and uj as well. The operator t* generates Zamolodchikov's 
integrals of motion mentioned above, while the asymptotic expansion of the function p sc with 
respect to its spectral parameter generates their vacuum expectation values. Another conjecture 
was that the asymptotic expansion of cj sc with respect to the spectral parameters describes the 
expectation values of descendants for CFT with central charge c = 1 — 6v 2 /(l — v) where v is 
related to the deformation parameter q = e mu . Equivalently, one takes the six- vertex model 
on a cylinder for some special boundary conditions and considers the scaling limit of the cor- 
responding normalized partition function with inserted quasi-local operator q 2aS ^0. On the 
other hand one computes a normalized three-point function in the CFT on an infinite cylinder 
with Virasoro descendants of some primary field <j) a of conformal dimension A a inserted at the 
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origin and two primary fields <f>± inserted at +00 and —00 respectively. The conjecture states 
that for any quasi-local operator one can find a corresponding descendant in such a way that 
the above partition function and the CFT three-point function are equal. In order to state this 
equality we used the fermionic basis for the quasi-local operators mentioned above. 

In [3] we established such a correspondence between the fermionic basis and the Virasoro 
generators only up to the level 6. It corresponds to bilinear combinations of the "Fourier- 
modes" j3j and 7^ with odd index j. For level 8 we have 5 Virasoro descendants l_ 2 , l_4l?_2> 
1?_4, l_6l-2> 1-8- As for the fermionic basis, we also find five linearly independent combinations 
of the creation operators /3* and 7*. One of them is quartic: fllPl-f^jl. Unfortunately, we could 
not uniquely fix the transformation matrix between the fermionic basis and the basis generated 
by the Virasoro descendants. We thought that it could be done if we would insert the simplest 
excitation corresponding to L—\<j>± at ±00 instead of the primary fields <f>± themselves. The 
difference between the "global" Virasoro generators L n and the "local" ones l n was explained 
in [3]. It is also discussed in the next section. 

The original motivation of the present paper was to fill this gap. It turned out, however, that 
it was not enough to consider only the simplest excitation. The above uncertainty still remained 
in this case. One has to take at least the next excitation, namely, that one corresponding to 
the descendants of the second level L 2 _ 1 cp± and L_2</>± in order to fix the unknown elements 
of the transformation matrix. All these elements are certain rational functions of the central 
charge c and the conformal dimension of the primary field A a They do not depend on the 
conformal dimensions A± of the primary fields <f>±. They are also independent of the choice 
of the excitation at ±00. This is rather strong condition. We are still unable to prove it for 
arbitrary excitation. We think that it is interesting to consider excitations also independently 
of the above concrete problem. Therefore we consider them in a more general setting and come 
to the solution of the level 8 problem in the very end. 

The paper is organized as follows. In Section [2] we remind the reader about the main results 
of the paper [JJ. Section [3] is devoted to the TBA approach for the excited states. We derive the 
equation for the auxiliary function O and define the function co sc in Section [U In Section [5] we 
discuss the relation to certain CFT three-point functions. We discuss the solution of the above 
level 8 problem in Section [6l In Appendix |A] we show several leading terms of the asymptotic 
expansion of the function $ defined in Section [3] and discuss its relation to the integrals of 
motion. In Appendices |B] and \C\ several leading terms for the asymptotic expansions of the 
functions F, x^, © are explicitly shown. 

2 Reminder of basic results of pEJ 

As was mentioned in the Introduction, in the paper [3J some specific connection between the 
conformal field theory (CFT) with the central charge c = 1 — 6^ 2 /(l — v) and the XXZ model 
with the deformation parameter q = e Ktu was established^! More precisely, the following relation 
was found to be valid with the left hand side containing the CFT data and the right hand side 
containing the lattice model data 



The left hand side means a normalized three-point function of the CFT defined on a cylinder of 
radius R parameterized by a complex variable z = x+iy with spacial coordinate x: —00 < x < 00 



(A_|P a ({U})<k(0)|A + ) 
(A_|0 a (O)|A+) 



lim 

n— >oo, a— >0, na=2nR 



Z K ' s {q 2aS ^0}. 



(2.1) 



^ook at the formulae (12.4) of [J or lpT8)) and (^9)) - (^Tll . 

2 Usually we take v in the region 1 < v < 1 called in [7] a "semi-classical domain" 
demands more accurate treatment. 



The region < v < \ 
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and coordinate in the Matsubara direction y: —irR < y < irR. The equivalence of the points 
x ± iriR is implied. At the origin z = some descendant field is inserted which is given by some 
polynomial P a ({l-k}) of Virasoro generators l_fc, k > acting on the primary field <j) a (z) with 
conformal dimension 

v 2 a(a — 2) , . 

= -jih^- (2 - 2) 

We called these Virasoro generators "local" in a sense that they are defined in vicinity of z = 
with the corresponding energy-momentum tensor 



T(z) = !> 



z- n -\ 



The bra- and ket-states |A + ) and (A_| are to be defined through two primary fields <f>± with 
conformal dimensions A± being inserted at x — > ±oo in such a way that L n \A + ) = 5 n fiA + \ A + ), 
n > when x = oo and (A_|L n = J nj oA_(A_|, n < at x = — oo. We called the Virasoro gene- 
rators L n "global" . They correspond to the expansion obtained via the conformal transformation 

z e~ z / R 



In [10] Alexander Zamolodchikov introduced the local integrals of motion which act on local 
operators as 



r dz 

i 2n _iO)H = / —h 2n (z)0(w) (n > 1), 



where the densities h,2n{z) are certain descendants of the identity operator /. An important 
property is that 

<A_|i 2n _ 1 (0(*))|A + ) = (/+_! - / 2n _ 1 )(A_|0(^)|A + ), (2.3) 

where I^n-i denote the vacuum eigenvalues of the local integrals of motion on the Verma module 
with conformal dimension A±. The Verma module is spanned by the elements 

iabi-i • • • iakp-il-aj • • • l-a, Wa(0)). (2.4) 

In case when A + = A_ the space is spanned by the even Virasoro generators {l-2n}n>l- 

In order to describe the right hand side of f)2. If) we need the fermionic basis constructed 
in [U[2] via certain creation operators. These creation operators called t*, b*, c* together with 
the annihilation operators called b, c act in the spac^l 

oo 

W (a) = W a ^ s , 

where W a _ S)S is the subspace of quasi-local operators of the spin s with the shifted a-parameter. 
They all are defined as formal power series of £ 2 — 1 and have the block structure 



t*(C) : W a _ SjS -> W 



a— s,s ; 



3 The problem of constructing the annihilation operator corresponding to the creation operator t* was discussed 
in the paper [13] but was not solved completely. 
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b*(C),c(C) : W a -.+i,s-l -> W 
c*(C),b(C) : W Q _ S _ 1)S+1 -> W 



a—s,s ? 



The operator t*(£) plays the role of a generating function of the commuting integrals of motion. 
In a sense it is bosonic. It commutes with all fermionic operators b(£), c(£) and b*(£), c*(£) 
which obey canonical anti-commutation relations 

[c(£),c*(C)] + = V(£/C,«), [b(0.b*(C)] + = -^(C/e,«) (2.5) 

with 

The annihilation operators b and c "kill" the lattice "primary field" q 2aS ^ 

k 



b«)(l 2 ° m ) = 0. c(O(9 2 ° s < 0) ) = o, S(k) = 



,4 



The space of states is generated via the multiple action of the creation operators t*(£), b*(£), 
c*(£) on the "primary field" q 2aS (°) . In this way one can obtain the fermionic basis. The 
completeness of this basis was proved in the paper [11] , 

In the right hand side of equation f)2. 1 1) we take the scaling limit of a normalized partition 
function of the six- vertex model on a cylinder with insertion of a quasi- local operator q 2aS ^ 



KSf 2aS(Q) , TVsTrM(y i i I - s) Ts,M^ S b^_ 1 ...b^ (^)o) 
Tr s Tr M (li- s) T s ,M q^ S • • • b^(g^W) N 



where the operators b^ • are defined through^ a singular part when £ — > and b* eg is a regular 
one 

c Q b*(c)po = ^r 2j b^-(x) + r a b; cg (C)W, x e w a _. +1>a _i 

i=o 

for some operator X of the spin s — 1, s > where C _Q b* eg (£)(X) vanishes at zero. The 
monodromy matrix Tg^ is defined via the universal i?-matrix of U q (sl2) on the tensor product 
of two evaluation representations where the first one fjs corresponds to the infinite lattice 
direction and the second one corresponds to the Matsubara direction 

rx -r^ 

oo n 

Ts,M = IJ Tj )M , Tj t M = Tj t M(X), 3j,m(C) = JJ Lj >m (() 
j=—oo m=l 

with the standard L-operator of the six vertex model 



Actually, later we discuss the Fateev-Dotsenko condition (|2.13[) fulfilled for the In this case 

the parameters a and s axe constrained. We will be mostly interested in the case when < a < 2. Then we 
need s < and another definition of the functional (|2.6p is necessary. It can be done through the replacement of 
operators b^j by c^-. This choice was taken in [15] with identification of the notation there c* creen _j = c^j. 
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The "screening operator" ^ carries spin —s. As was discussed in [3], the functional (|2.6p 
does not depend on the concrete choice of the screening operator under rather mild conditions. 
Due to common wisdom, in case of an infinite lattice, one can change boundary conditions and, 
instead of taking the traces in the right hand side of (|2.6I) . insert two one-dimensional projectors 
and \k + a — s, s){k + a — s, s\ at the boundary, where \n) is eigenvector of the transfer 
matrix Tm(C> k ) = T^j (Tj t MQ Kai ! ) with maximal eigenvalue T(£, k) in the zero spin sector and 
where the eigenvector \k + a — s, s) corresponds to the maximal eigenvalue T(£, k + a — s, s) of 
the transfer matrix Tm(Cj k + a ~ s ) m the sector with spin s. The twist parameter k plays the 
role of the magnetic field. So, we can perform the following substitution 

{q 1 (K + a-s,s\T SMq ^s h ^ s _ i ... h ^ q 2 a s(o ))lK) > 

which does not affect the answer for the case of infinite lattice if 

( K I^M \ K + a ~ S > S ) ^ 0- 
The theorem proved by Jimbo, Miwa and Smirnov [3] claim!] that 

Z«> s {t*(Q(X)} = 2p(C|K, « + a, S )^ S {X}, (2.7) 



Z K > s {b*(C)(X)} = _L jf w ( C ^|K,a,,)^{c(0(X)}^, 
Z^{c*(C)(X)} = -^jfa;(^,C| K ,a, S )Z^{b(OW}^, 



where the contour T goes around all the singularities of the integrand except £ 2 = £ 2 . The direct 
consequence of the above theorem and the anti-commutation relations (|2.5p is the determinant 
formula 

Z^fK?) • • • t%°)b*(d + ) • • • b*(C+)c*(C r -) • • • c*(Ci)(Q 2am ) } (2-8) 
v 

1=1 

The functions p and a; are completely defined by the Matsubara data. The function p is the 
ratio of two eigenvalues of the transfer matrix 

/>, . \ T((,k + a - s,s) 
p{QK + a-s,s) = — . 2.9) 

We will come to the definition of the function oj in Section HJ in more general case of presence 
of the excited states. 

The scaling limit in the Matsubara direction means 

n -> oo, a — > 0, na = 2ttR, (2-10) 

where the radius of the cylinder R is fixed. Simultaneously one should rescale the spectral 
parameter 

£ = Act", a = Ca (2.11) 



Actually, in [3] the statement was proved for the case s = but as was discussed in [4], this statement can 
be proved for s =fc also. 
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with some fine-tuning constant C. One of the most important points of [1] was to define the 
scaling limits of p and u> 

P sc (X\k,k')= lim p(Xa u \K,a,s), (2.12) 

scaling 

u sc (X, p\k, k' ,a) = — lim u}(Xa u , pa u \n, a, s) 

4 scaling 

where n' is defined through an analogue of the Dotsenko-Fateev condition |16j 

k! = k -\- ct -\- 2————S. (2.13) 
The continuum limit can be taken in both directions of the cylinder. 

The first conjecture proposed in [3] was that the creation operators are well-defined in the 
scaling limit for the space direction when ja = x is finite 

r*(A) = limV(A^), P*(X) = limib-CAo"), J*(X) = IimJc*(Aa") 

a— >0 Z a— >0 Z a— >0 Z 

and for the "primary field" 
$ a (0) = lim q 2aS (°). 

a— >0 

Asymptotic expansions at A — )• oo look 

CO 

log(T*(A))~ XX-iA- 2 ^, (2.14) 

oo 1 oo 

77ie nexi conjecture based on the bosonisation arguments was that the scaling limit of the 
space W a _ SiS belongs to the direct product of two Verma modules 

Scaling limit (W a _ s>s ) C V |2 i-< <g> V_ Q 

and the operators t*(A), /3*(A), 7* (A) act non-trivially only on the first chirality component 
and do not touch the second one 





-1 : 






/% 


-1 : 


V a+2i^(s-l) 




72,- 


-1 : 


V a+2i^(s+l) 





Using the results of the paper [7], we get the asymptotic expansion for A — > +00 

00 

logp sc (A| K , K ') ^X-^C j (l+_ 1 - (2.15) 
i=i 

Here the integrals of motion l^j—x are ^ ne same as m (|2.3p . They correspond to the "right" and 
the "left" vacuum and depend on k and k' respectively: 
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We can identify 

For the function d sc we have 

oo 

UJ SC (X, fl\K, k' , a) ~ \J P sc (X\k, K , )y / p sc (fi\K, K 1 ) ^ A ~ [I J ~LUij(n, k' , a) (2.16) 

*>J=1 

when A 2 ,/i 2 — > +00. 

The scaling limit of (|2.8p is proportional to 

Z^'{r*(\\) ■ ■ ■ r*(A£)/3*(A+) ■ ■ ■ /3*(A+) 7 *(A-) • • • 7*(AD(*«(0)) } (2.17) 
p 

= Y[ P sc (\°\k,k') xdet ( W sc (A+,AT| K , K ',a)) M=w . 
i=l 

If we substitute the expansion (|2.14p into the left hand side of (J2.17P and the expansions (|2.15p , 
(|2.16p into the right hand side of (|2.17p we can compare the coefficients standing with powers 
of the spectral parameters and express the functional Z^ K of any monomial constructed from 
the modes t\^_ x , fi*2j-ii 72j-i through the integrals of motion hn-i, coefficients C n and u)ij. 
In [1] we argued that the eigenvalue T(£, k + a — s,s) in the scaling limit (|2.10p becomes equal 
to T(£, k'). This means that, if we choose a and the spin s in such a way that k' = k, then 
using (j2.9[) and (|2.12p . we get p sc (\\k, k') = 1. This is an important technical point. In this 
case we were able to apply the Wiener-Hopf technique and obtain the coefficients as power 
series in k~ 1 where k — > 00. 

On the other hand, one can evaluate the left hand side of (|2.1|) using the operator product 
expansion (OPE). In order to compare with the result obtained by the above lattice method, 
one needs to identify the parameters k, k', a with the CFT data. In fact, we already identified a 
by taking the formula (|2.2p for the dimension A a of the primary field <p a . The next step was to 
take A + = A K+ i and A_ = A_«/+i. 

The most important conjecture of [4] was that for the CFT with central charge c = 1 — 
it is always possible to find one-to-one correspondence between a polynomial P a {{\-k}) in the 
left hand side of (12. lj) and certain combinations of /Sjy-i an d Tsy-i- 

It is convenient to introduce 

Plm-l = £> 2m-l(")/32m-l> 72m-l = Ajm-l(2 - a)72m-!> 

where 

1 r(f + i(2m-l)) 



, , 1 , , 2m-l , , 2m-l 

£>2m-l(a) = -75=r(l/) ~(1 " f 



(ro-l)! r (§ + (^)( 2m _i) 



, 2 1 2i/ 

together with even and odd bilinear combinations 

1 



'2771-1,271- 



-1 - (m + n- 1)- (P 2m _i72n-1 +^2n-l72m.-lJ ' (/.IS) 



^odd , _ in 1 /aCFTi. CFT* /oCFT* CFT*\ 

^2771-1,271-1 — "a l"i-T — -LJ - ^P2n-1 72777.-1 ~~ P2m-l72n-l ) i 



where 



v{y-2) 



d a = ^_ i J (a-l) = - V(25 - c)(24A a + 1 - c). (2.19) 
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The Verma module has a basis consisting of the vectors (|2.4p . Conjecturally the same space is 
spanned by: 

l2fei-l • • • i2fc p -102mi-l,2m-l • • • 02^-l,2n r -102mi-l,2fii-l < ^2^-l,2n fI -l(^Q!)- (2.20) 

Since 

[ IciVJ = (2j - l)r^_ l5 [ lo,/35i_i7^--i] = (2« + 2j - 2)/%_ l7 £ i _i, 

the descendants of the form (|2.4p and those created by the fermions of the form (|2.20p must be 
finite linear combinations of each other if the corresponding levels coincide 

q r f 

Yl 2l i = Yl 2 ( mk + nk-l) + J2 2 ( m £ + - !)■ 
j=l fe=l fc=l 

As was discussed above, one can choose the parameters a and s in such a way that n' = k and 
A + = A_. With this choice we factor out the integrals of motion. The quotient space of the 
Verma module modulo the action of the integrals of motion is spanned by the vectors of the 
form 



Uh ■■■!-«, (0«(O)). 

All coefficients of the polynomial P a ({l_/ : }) are independent of k. We were able to identify 
the above basis vectors up to the level 6. The system of equations is overdetermined but 
nevertheless it has a solution: 

0S7^U, Cf = l 2 _2 + ^^l-4, <3 d = ^-4, (2.21) 

-even „ ,3 C + 2 - 20A + 2 C A 



3(A + 2) 

5600 + 428c - 6c 2 + 2352A - 300cA + 12c 2 A + 896A 2 - 32cA 2 



60(A + 2) 

odd^ 2A 56-52A-2c + 4cA 

i.» = aT^ 1 - 41 - 2 + sTaT^ u ' 



l 



-even ^,3 , 6 + 3c-76A + 4cA 

-6544 + 498c - 5c 2 + 2152A - 314cA + 10c 2 A - 448A 2 + 16cA 2 
+ 60(A + 2) ~ 6 ' 

where for simplicity we took A = A a . We hope it will not cause some confusion by mixing 
this A with the anisotropy parameter of the XXZ model. In the above formula (|2.2ip we imply 
only a weak equivalence "=" between the left hand side and the right hand side. The sign = 
means that the left and right hand sides being substituted into the corresponding expectation 
value give the same result by acting on $ a (0) and (j) a (0) respectively. In other words A = B if 
and only if 

KKt „/ , xxX <A_| J B(0 Q! (O))|A + > 



(A_|0«(O)|A. 



with the functional Z^ K defined in (|2.17|) . In the next sections we will generalize the scheme 
of [1] to the case of excited states. 
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3 TBA for excited states 

In fact, the Jimbo, Miwa, Smirnov theorem (|2.7|) is valid not only for eigenvectors \k), \k+o.— s, s) 
associated with maximal eigenvalues in the corresponding sectors, but also for eigenvectors 
corresponding to excited states. Here we consider only special excited states, namely, the so- 
called particle- hole excitations [17]. Some motivation for this is as follows. In [7J Bazhanov, 
Lukyanov and Zamolodchikov formulated several assumptions about the analytical properties 
of the eigenvalues of the Q-operator with respect to the square of the spectral parameter £ 2 , in 
particular, that its zeroes in the complex £ 2 -plane are either real or occur in complex conjugated 
pairs, that an infinite number of zeroes are real and positive and that there may be only a finite 
number of complex or real negative zeroes. An important assumption is that the real zeroes 
accumulate towards +oo in the variable Q 2 . Further, in [TH] the authors proposed a conjecture 
that for the asymptotic analysis, when the parameter n becomes large, only real positive zeroes 
corresponding either to the vacuum or to the excited states are important. The question why 
real negative zeroes or complex zeroes are not important for the asymptotic analysis seems hard. 
Usually one uses experience coming from the numerical study and also the arguments based on 
the analysis of a small vicinity of the free-fermion point v = 1/2. Counting arguments play 
important role here as well. But we do not know any rigorous proof of this statement. Any 
further discussion of this subtle question would lead us beyond the scope of this paper. 

Let us start with the case of a lattice with an even finite number of sites n in the Mat- 
subara direction. The Bethe ansatz equations (BAE) are usually deduced from the Baxter's 
TQ-relation [9] for the eigenvalues T(£, k, s), Q(£, s) of the transfer matrix Tm(C> k ) an d the 
Q-operator Qm(C> k ) defined by the formula (3.3) of [J| 

T(C, K, s)Q(C, K, s) = d(()Q(q(, k, s) + a(C)Q(q- l C, K, s) , Q((, k, s) = C K+S A((, k, s) 

with polynomial dependence of the functions T and A on the spectral parameter £ in every spin 
sector s, < s < ^ and 

«(0 = (i-?c 2 ) n , d(o = (i- q - 1 c 2 ) n - 

In the TBA approach one introduces auxiliary function 

rf(C)Q(gC,«,s) 

which satisfies the BAE 

n 

a(£j,K,s) = -1, j = l,...,--s (3.1) 

with § — s zeros £j of Q(Ci k, s) called Bethe roots. The BAE in the logarithmic form look 

log a(£j, k, s) = nirrij, (3-2) 

where rrij are pairwise non-coinciding odd integers. The ground state corresponds to s = and 
the choice rrij = 2j — 1. The following picture describes this situation 

•••oooo»»»»»--- 

-5-3-1 1 3 5 

where the black circles correspond to "particles" and have positive odd coordinates rrij while the 
white circles correspond to "holes" and have negative coordinates. The particle-hole excitations 
can be got when some of "particles" are moved to positions with negative coordinates. Let 
us denote /(+> fc ) an ordered subset of positive coordinates l[ + ' k) < ■ ■ ■ < < n which 
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correspond to the positions of created holes and 1^ correspond to negative coordinates of the 
moved particles^ —ii < ■ ■ ■ < —l[ ' k ^: 



■5-3-11 3 5 ■ 



We will denote the corresponding vector and co-vector \k; l( + ' k \ /( and (k; /( + > fc ), /( > fe )|. 
Let r = 1, . . . , k be the Bethe roots corresponding to the particles moved into the positions 
with negative coordinates —Ir and correspond to the holes with positive positions Ir + ' k \ 
Following the papers [5j [T9| [20] , we can rewrite the B AE f|3. 1 1) in form of the non- linear 
integral equation 

log o(C, «, s) = -2tt^(k - a) + log f " / log (1 + a(£, «, s)) , (3.3) 

where the contour 7( S) fc) goes around all the Bethe roots £j including the moved ones £j r = £~ 
in the clockwise direction and the kernel 

K((, a) = ^A C ^(C, a), K(Q = K((, 0), A C /(C) = f(qQ ~ f(q~\)- 

For the case of the ground state s = 0, k = we have 7(o,o) = 7 with the contour 7 used in 
Section 3 of [I]. Let us for simplicity stay with the case s = 0. One can rewrite (|3.3p replacing 
the contour 7(o,fc) by 7 and taking into account the contribution of the residues from the moved 
Bethe roots £~ and holes £+. There are also other equations coming from (|3,2p . Altogether we 
have the following set of equations 



k 

log o(C, k) = -2niu K + log + £ {9((/Zr) - <7(C/£T)) 

logo(^,«) = T7ri/ r (± ^, r = !,...,£, (3.4) 



-co; 



«-(c/oiog(i + o(e,«)) ^ 



1 2/-2 

where 5(C) = log jf^g • 

Now we can consider the scaling limit and generalize the analysis of Section 10 of [1]. Let 
us start with a small remark. So far, we considered only particles which "jump" to the left. In 
principle, for the case of finite n we should also consider particles which "jump" to the right. 
But in the scaling limit when n — > 00 it is implied that the right tail of the Bethe ansatz phases 
in the ground state is infinite. So, the jumps to the right are irrelevant in the scaling limit. 
Therefore we will not discuss them here. Let us introduce the functions 

T sc (X,k)= lim T(AaV), 

n— >oo, a— >0, 2irR=na 

Q sc (A,k)= lim a UK Q(Xa u ,K) 

n— >oo, a— >0, 2ivR=ria 

and taking into account that for 1/2 < v < 1 the ratio a(£)/d(Q — > 1 in the scaling limit with ( 
related to A as in (|2.1ip . we get the scaling limit of the auxiliary function a as 

a (A, k) - 



Q sc (Xq-\KY 



The case of s > can be treated similarly to the case s — 0. The only difference is that for the finite n case 
one has less Bethe roots. 
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We would like to study the asymptotic behavior of o sc for A 2 — > oo and k — > oo in such a way 
that the variable 

is kept fixed with 

c{v)=T{v)- 2 e s [^f\ 5 = -u\ogu-{l-u)\og{l-u). (3.5) 
Then the function 

F(t, k) := loga sc (A,K) 
satisfies the equation 

f 00 du k 

F(t, k) - / K(t/u)F(u, k)— = -2-nivn + V(t?(t/ti") - g(t/t~j) (3.6) 

- / log (l + e F ^) — + / tf(t/u) log (1 + e - F ^)-, 

Ji u Ji u 



where e is a small positive number and 
t± = c („)-i 



v ' {no) 2 "' 
With a slight abuse of notation we will write 

1 1 /^2 + l tq~ 2 + l \ tq 2 -l 

Introducing the resolvent R(t,u) 
f°° dv 

R(t, u) - / —R(t, v)K(v/u) = Kit/u) (t, u>l) 
Ji v 

and 

G(t;u,v) = G(t,u) - G(t,v), 

f°° dv 

G(t,u) = ((I + R)g)(t,u)=g(t/u)+ R(t,v)g(v/u) — , (3.7) 

Ji v 

one can rewrite the equation (|3.6p as follows 

k 

F(t,K) = nF (t) + ^G(t;i+ t") 

r=l 

- R(t,u)\og(l + e F ^)-+ R{t , u )log{l + e- F ^)-, (3.8) 

Ji u Ji u 

where if depend on k and can be defined from the equations 

F(tf,K) = Tm4 ± ' k \ r = l,...,k (3.9) 
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and Fq is the same as in [3]. It satisfies the integral equation 

((/ - K)F )(t) = -2mv, (3.10) 
where as in ()3.7p . the contraction is defined on the interval [1, oo) 

Kf(t)= [°° K(t/u)f(u) — . 
Ji u 

It follows from the WKB technique [7] that the asymptotic behavior at t — > oo 

F (t) = const • i 277 + 0(t"277). 
Then the equation (|3.1U|) can be uniquely solved by the Wiener-Hopf factorization technique 

F (t)= [ dlt u S{l)—^- (t>l), 

JR-J--iO 1(1 + or:) 



where 



/ = — ; 1 (3-11) 
2^/2(1^) 

and the function 

= T(l + (l-u)ik)T(l/2 + tuk) piSk 
T(l + ik)y/2ir(l-v) 



with 5 defined in (|3.5p satisfies factorization condition 

i - = s(k^s(- k )-\ k( k ) = sinh(2 ;- 1)7rfc , 

smh7rK 

where K stands for the Mellin transform of the kernel K 

f'°° rlt 1 f'°° 

K{k) = / K{t)t~ ik — , K{t) = — / K(k)t ik dk. 

J0 t 27r J-oo 

Below we also use the asymptotic expansion of S(k) at k — > oo 

oo 

S(k) ~l + J2Sj{ik)- j . 

3=1 

For example, the two leading terms look 

_ (l + u)(2u-l) _ {l + vf{2v-lf 

bl ~ 24(1 -!/> ' 2 " 1152(1 - V )W ' [6A2) 

The same method leads to 



-oo 27r J-oo 2?r / + m - zO 



/ ™ 1«u™S(l)S(m)K(l)K (m) — -. (3.13) 

.oo 27r J.oo 27T t + m - ^0 
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One can solve the equation ()3.8[) iteratively using the asymptotic expansion 



F(i,re) = ]T, 



-n+l 



F n it). 



n=0 



An important difference of this expansion with the formula (10.3) of [1] is that here also even 
degrees of re may appear. Since 

-^- u ^(t/u) = K(t/u), 

we can get the following equation by differentiating (|3.7p 

-^-u^G{t/u) = ((/ + R)K)(t, u) = R(t, u) 
2m ou 



and then using (13.13P 



G(t,u) 



-2iri 



dl 



2tt 



dm 



£ / ^1 fu™S(l)K{l)S{m 



—i 



—i 



2tt 



I + m — iO m + iO' 



where we used the regularization m + iO in such a way that G(t, u) — > for u — > oo but in fact, 
this does not matter because G(t, u) enters into the equation (|3.8jl only through the difference 
G{t,u) -G(t,v). 

One introduces the function ^(1,k) which has an asymptotic expansion 



if 



n=0 



(3.14) 



related to the function F(t,n) via 

/oo 
dl t u S{l)K(l)(V(l, re) - k9 q (1)). (3.15) 
-oo 

It is convenient to introduce 
P ■= /«, 

where / is defined in (|3.1ip and consider the asymptotic expansion with respect to p instead 
of re. 

Using (|3.15j) . one can rewrite (|3.6p in the following equivalent form 



r=l 



P ./0 



^^,e^)log(l + e F(p) ( e " /p ^) 



+ J g/J(l >e — /*)Iog(l + e- 



) ■ 



(3.16) 



where e is a small positive number and by definition 

F«(t,p) :=F(t,re), *W(|,p) :=*(/,«), *«(Z) := 1* (0 = -—J-^ 



and then 

R(l,e ix / p ) = Tea h 



e -hx/p 

l + h 



■S(h) 



G(l;y,z) = -res h 



-hy/p 



-hz/p 



h{l + h) 



-S(h) 
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where vesh is the coefficient at h 1 in the expansion at h = oo. Taking 

F(p)(e ix / p ,p) = -2vr (x - F(x,p)) , 
we obtain the Taylor series at x = 

F(x,p) =x + res h [e- hx / p S(h)i¥P\h,p)]. (3.17) 

As the asymptotic series with respect to p at p — > oo, it starts with p~ l . The 2k parameters tf- 
participating in (|3.9p and f|3. 16|) are functions of k or equivalently of p. We take them in the 
following form 

t±(«) = x ±(p) = f>±p- J '. (3.18) 

j=0 

Then we get from (|3.16|) 

fc 

*V(l,p) = £G(Z;:r+(p),*-(p)) - —^— +H(l,p), (3.19) 

where 

X E^(^) n (-£j I l°g(l + e" 2 ^). (3.20) 

n=0 ^ even 

Here {f(x, J|)} C vcn = \{f{ x i -§^) + f(—x, — g|))- The parameters x^(p) are to be determined 
from the condition 

F(xt(p),p)=xt( P )T^ ± ' k) (3-21) 
which is 0(p^ 1 ) if 



.r 



, ± o = ±^4 ± ' fc) - (3-22) 



r,0 2 

So, we come to the iterative scheme which allows us to compute F^ p \t^p) to any order of p~ l . 
One can see that the following expansions are consistent with the above equations: 



F(x,p) =Y^F j {x)p- j -\ (3.23) 

3=0 

oo 

H(l,p) = J2H j (l)p~ j - 1 , (3.24) 



j=0 

oo 



G(l;y,z) = Y J Gj(l;y,z)p- J - 1 (3-25) 

3=0 

and in agreement with (I3.14|) 

oo 

^\l,p) = Y^^>f{l)p~ j+1 - (3.26) 

3=0 
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Let us explain the very first iteration step. We first compute Hq in the expansion (|3.24|) 
using (|3.2U|) . the formula 



oo 



^ x m (9\ n log(1 + e -2n X) = m]{1 _ 2 -m-i + n, C(m n + 2) 



2vr V dx J K / -v y ( 2 7r) m - n+2 

and the fact that only the term with n = in the sum at the right hand side of (|3.2U|) contributes 
because F is of order p" 1 . So, we easily come to the result 

*mo = -i- 

Then we substitute it into (|3.19j) and calculate the leading order of the function G taking into 
account the condition (|3.22j) . As a result we obtain a few leading orders of ^>&>; 

^p) = -^fxy + (E ^(4 + ' fc) + 4~' k) ) - ^ p- 1 + o(p-). 

This we substitute into the formula f)3. lTj) and get Fq from (|3.23p 



where Si is given by (|3.12p . Then we take the equation (|3.2ip up to the order p 1 and easily 
solve it 



/ A±,k) \ 2 f k 



One can do the second iteration by repeating this procedure. In Appendices |A] and [B] we show 
the result of such a calculation for few next orders with respect to p^ 1 . 



4 The function uj sc 

Also we need to generalize the expressions (11.5), (11-6) of [1] for the function lo sc to the case of 
excited states. Still we take the condition k = k' for which we can choose the excited state for 
the spin sector and for the sector with spin s in such a way that p(X\n, k') = 1. The reasoning 
here is quite similar to that one described in Section 4 of [4]. Actually, one can start with the 
similar expression to (11.1) of [1] but with a generalized dressed resolvent -Rdrcss 

W SC (A, fl\K, K, a) = (/i eft * fc / ri ght + /left *fc -Rdress *k /right) (A, fi) + W (A, fJ,\a), (4.1) 

wher^l 

/left (A, M> «) = 2~^^a^o(A///,Q!), /right (A, M> «) = S~ij}o(X/fi, a), 
<5a/( A ) = /(?A) - /(A), wo(A,/i|a) = S^A^^X/fra), 
^o(A, a) = A A /(A) = /(gA) - /(g^A) 

7 In comparison with (11.1) of [4] we take instead of the function ip defined in (|2.5|) the function ipa. The result 
does not change if we also change the kernel K a — > K a , Q as it is done in (|4.2[) . 
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with 



A A Vo(A,a = -VP / — rr^r - T' 



where the principal value is taken for the pole p? = 1. The contraction *fc means 

dA 2 



f *k 9= f{X)g(X)dm(X), dm{\) = — 

with the contour 7(o,fe) which corresponds to 7(o,fc) from the equation (I3.3h but taken for variab- 
le A 2 instead of £ 2 = X 2 a 2u . We remind the reader that the contour 7(o,fe) was taken around 
all the Bethe roots in the clockwise direction in case of the excited state with k particles and k 
holes for the zero spin sector. 

The resolvent i?dress fulfills the integral equation 

-Rdrcss - -Rdrcss *fc K a ,Q = K a>0 , K a>0 (X) = A A ^ (A,a). (4.2) 



Applying a similar trick which we used to derive (|3.4j) . namely, deforming the contour 
7(o,fc) — ^ 7(o,o) and taking into account additional terms coming from the residues corresponding 
to particles and holes, we can obtain 

D u \ ou ^ o • f R{t,t+,a)R d Tess(tt, u ) R(t,t~ ,a)R diess (t~ ,u) 
Rdrm^(t, u) — Kit, u, a) = Zin > —r 



r=l v 

R(t, v; a)R dress (v, u) dv f e lE °° R(t, v; a)R dress (v, u) dv 



1 + e - F ^' K ) v Ji 1 + 



(4.3) 



F'(t,K) :=t^F(t,n), 



where as in [4], we introduced the "undressed" resolvent R(t,u,a) which satisfies the equation 
f°° dv 

R(t,u,a) — / — R(v, u, a)Ko(t/v, a) = Ko(t/u, a) (4-4) 



l v 



with the kernel 



Ko{t,a) = — 



2ttz \tq 2 - I tq~ 2 -11 

corresponding to the above kernel K a $. The solution of (|4.4[) again can be got by the Wiener- 
Hopf factorization technique 

R(t,u,a) = Ko(t/u,a) 

/°° dl f°° dm —7 
— / — t il u im S(l,a)S(m,2 - a)K(l,a)K(m,2 - a)- 
.oo 2vr J_ OQ 2tt 1 + m-iO 



with the Mellin-transform 

sinhvr((2^ - l)Jfe - f ) 



K(k,a) 



sinh7r(A; + f ) 



corresponding to the kernel Ko(t,a) and 

_ r(i + ( i-^ f )r(i + ^) 
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1-K(k,a) = 5(A:,a)- 1 5(-A:,2-a)- 1 . 
We assume that 

< a < 2. 
Now we take the ansatz for Rd ress (t,u) 

Rdrcss(t,u) = K (t/u,a) (4.5) 

/ t il u im S(l,a)S(m,2 - a)K(l,a)K(m,2 - a)@(l,m\p,a) 

-oo 27r J-oo 27r 

with the asymptotic expansion^] at p — > oo 



@(l,m\p,a) ~^2@ n (l,m\a)p 11 , @o(l,m\a) 



n=0 



I + m 



(4.6) 



and substitute it into the equation (|4.3j) . As a result we get the equation which allows us to 
calculate all 6 n by iterations 



m\p, a) — Qq(1, m\a) = - resj/res^ 

/ k k 

(l'+m')xt(p)/p („\ _-\ _ -[■\iA\^ t ,-( l '+ m ') x r(p)/p/(p' 



J2 e -(i'+m')x r (p)/p/(f'(x+( P ),p) -l) + iY j e 

r=l r=l 
n=0 n ' ^° ^ 



X — « 2^ e 

V r=l 

oo 



S(Z', 2 - a)S(m', a)Q{m', m\p, a) /{I + V) 

/(F'(x;(p),p)-1) 

(4.7) 



dx 



1 + e 2 ™ 



where the odd part {f(x, ^)}odd = \{f{%, ^) - and 
d - 

F (x,p) := g-F(x,p). 

Performing iterations implies that enough many orders in the expansion of F(x,p) and x^(p) 
with respect to p~ l were obtained by means of the iteration scheme described in the previous 
section. For a few leading terms in the expansion for Q(l,m\p, a) we get 



m\p, a) 



l + m \ 24i> 2v 



r=l 



\-iu(l + m) - i + A a ^j p~ 2 + 0(p~ 3 ) , 



where A a is given by (|2.2p . Some other terms of this expansion will be shown in Appendix ICl 
for the case h = 1. 

One can derive relation like (11.5) of [1] using the form (|4.ip . integral equations (|4.2p . (|4.3p 
for the dressed resolvent -Rdress and the ansatz (14.5P 



w sc (A, /ti|/c, k, a) 



(4- 



8 Here Q n (l, m\a) are different from those introduced in [4] since we should take into account the contribution 
of terms with odd degrees n also. 
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dZ / dm 5(Z, a)S(m, 2 — a)0(Z + iO, m\p, a) 



~ 2iri J—oo J —oo V ' ' V ' y V ' 'V* 2 ^)/ V« 2 ^) 

where we returned to the variables A, and 

r(-<fe + f)r(i + ti/A;) 



5(/c, a) 



r(-*(i - v)k + § ) V§7r(i - v)^- 



■a)/2 ' 



The asymptotic expansion at A, fx — >■ oo can be obtained by computing the residues of the 
functions S(l, a) and S(m, 2 — a) 

oo ^ 

E2r— 1 2s— 1 

— ■ -D 2 r-i(a)D 2s -i(2 - a) A ~/z ~fi 2 r-i,2 S -i(p, a),(4.9) 
r + s — 1 

r,s=l 



where 



1 , N 2n-i, 2n-i 1 r(f + ±(2n- 1)) 

vV (n-l)! r ^ + iigd(2n-l)J 

and 

\ 2r+2s-2 



The relations (|4.8f> — (|4.11)) look the same as (11.5)— (11.7) of [3J. However there is an important 
difference. It was pointed out in [4] that the expansion coefficients n ( ^^D ^ ' '^iL ^ P' Q ) 
satisfy the so-called vanishing property i.e. for given r and s they vanish starting from 2n = r + s. 
It is equivalent to the fact that the function G is proportional to a polynomial with respect 
to p. This is true only for the case of ground state k = and also for the case k = 1 with 
j(+,i) = = i. For both cases all the coefficients n with odd n vanish. We will see that 

in both cases the space of the CFT-descendants is one-dimensional. 



5 Relation to the CFT 

Here we would like to study generalization of the relation (|2.ip between the lattice six vertex 
model and the CFT to the case of excited states. As discussed in Section [2J for the case 
of the ground state we inserted at ±oo-points of the cylinder the two primary fields <f>± with 
dimensions A±. We also identified A + = A K +i, A_ = A_ K / + i and A + = A_ for the case k — k . 
The corresponding states were denoted |A_|_) and (A_| respectively. For the six vertex model 
we introduced in Section [3] the states \k; I^ + ' k \ /(~> fc )) marked by the two ordered sets 
and of k odd, positive, non-coinciding numbers. In case of the spin s sector we denote 

such a state \k, s; l( + ' k \ In the scaling limit we can identify these states with descendants 

of the primary fields <f>± at level N = | ^ (Ir + ' k ^ + Ir ' k ^) 

r=l 

|«; JC-.*)) -+ scal IA+j/C+^./C-.*)) = £ AC' k ^~' h)) L-ni---L-n m |A + ), 

n-x>"->n m >l 
n 1 -\ \-n m =N 

( K + a -s,s;I^ k \l^\ ^ scal <A„;/(+' fc ),/(-> fc )| 

= E A(£: k ^-' k) \A^L nm ...L ni . (5.1) 

Th\>--->nm>l 
ni-i \-nm— N 
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We normalize 

(A_;j( + ' fe ),j(- fe )|0 Q (O)|A + ;j( + ' fe ),/^) = 1 (5.2) 

but we do not demand orthogonality of the states with different sets. In other words, the "scalar 
product" (A_;/( + ' fc ),/(-' fc )|0 a (O)|A + ;/ /(+ ' fc) ,/ /(_ ' fe) ) is not necessarily zero. So, instead of ([2TTD 
wc take 



/ ; / ^ni,...^™ n' ...,n' , 



nj^^-'-^nm^l n ' >...> n ' >l 

ni+-+n ro =JV 1 ™' 

n' H hn ,=JV 

— ni ' ' ' L—rim |A+) (5.3) 

(k + a - s, s; /(-' fc )|Ts,Mg 2K5 b^ s _ x • • • (q 2aS ^O) \k; /(+-*), /(-*)) 

= lim ■ ■ 

™* <« + « - s; ^ (+ ' fc) ^ ( -' fc) |r S ,Mg 2K5 bJ 0jS „ 1 • • • b^ )0 (g2aS(0)) | K . /(+,*),/(-,*)) • 

na=27r_R 

In this relation the polynomial P Q ({l_j}) does not depend on k and the choice of the excitation 
i.e. it is independent of k and of the both sets I^' k ^ . As was pointed out above, the coefficients of 
this polynomial are rational functions of the conformal charge c and the conformal dimension A a 
of the primary field (j) a only. If the operator = $?L\ ■ ■ ■ P^Z-ilifJ-i • • • I2I-1 then the 

n 

level of descendants participating in P a ({l_j}) is M = 2^(ij + ji — 1). On the other hand, 

1=1 

the coefficients A nii ___ >rim , A nij ... jnm are independent of the choice of the operator 0. So, our 
strategy is to take for a given level M all linear independent operators represented in terms of 
the fermionic basis modulo integrals of motion and as many different excitations as necessary in 
order to fix the corresponding polynomials P a ({l_j}) and the coefficients^ A nir ,, j7lm , A nir „ >nm . 
All other relations for this level partially determine further coefficients A, A and the rest of 
the equations fulfills automatically. In [4j we were able to fix the polynomials P a ({l-j}) up to 
M = 6 using only the ground state data when k = 0. It means that for M < 6 the relation (j5.3[) 
fulfills automatically for any excitation i.e. for any k and any two sets I^^ k \ We checked this 
for the case k = 1 with = = 1 and = 1, = 3, = 3, = 1. 

Starting with M = 8 the situation changes. We do not have enough equations in order to fix 
the polynomials P a ({l_j}) if we restrict ourselves with the case of the ground state. We need 
additional equations involving excitations. In the next section we will consider in detail the case 
M = 8. 

Before we go further let us make one remark. It is interesting to note that the above fermionic 
basis operators are marked in exactly the same way as the particle-hole excitations, namely, by 
two ordered sets of n odd, positive integers j(+' n ) = {2i\ — 1 < • • • < 2i n — 1} and /(~ ,n ) = 
{2ji — 1 < • • • < 2j n — 1}. So, we can denote such an operator 

n - aCFT* ^CFT* CFT* -,CFT* /r a \ 



In case we take this operator in the right hand side of (|5.3p we will use the following shorthand 
notation 

r.h.s. of Q := (/ (+ ' fe) ,/ ( -' fe) |0 /C+ ,„) jjC -, n) |/( + ' fe ),/(- fe )}. 

Now let us describe in general how we compute both sides of the relation (|5.3p . We start 
with the right hand side of (I5.3P determined by the lattice data. We pointed out the fact that 



9 In fact, not all coefficients A, A can be fixed but rather some of their products. We checked for one particular 

case that one can fix them completely taking A ni: ... t „ m = A ni „ m but a'priori it is not quite clear to us why it 

should be so. 
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the theorem by Jimbo, Miwa, Smirnov works for the case of excited states also with the same 
fermionic operators. Only two transcendental functions p and oj are sensible to the changes that 
happen in the Matsubara direction. It means that we still can apply the Wick theorem and after 
taking the scaling limit come to the same determinant formula (|2.17p but with a new functional 
determined by the right hand side of (|5.3[) instead of Z^ K . With our choice of parameters a, s, 
k, k' the function p sc is still 1 and the function u sc is now determined through the asymptotic 
expansion (|4.9j) . So, we come to 



(/(+,*), j(-,fc)| 0j(+in)j/( _ in) |j(+,fc) j/ (-,fc)) = det ( 



^2v-l,2j r ,-l(P,a) , 

(5.5) 



i r + j r > - 1 



r,r — l,...,n 



where the function is defined by (|4,lip . 

Now let us proceed to the left hand side of (|5.3p which involves the CFT data. For simplicity 
let us put the radius of the cylinder R = 1. We need to calculate 

<A_|L n / • • • L n >P a ({\^})(j) a {Q)L^ ni ■ ■ ■ L_ n JA+) 
with n\> ■ ■ • > n m > 1, n'-y > • ■ ■ > n' , > 1 and 



m 

j=l 3=1 



or picking out some monomial with respect to the local Virasoro generators we need the 
following value 

Ki , • • • ) n' m , ; ai, . . . , a d ; m, . . . , n m ) 

■= {A-\L n > ■■■ L n > l-2a d ■ ■ ■ l-2oi<^a(0)i- m ■ ■ ■ L_ nm \A + ) (5.6) 

m 1 

with d positive integers 1 < a\ < ■ ■ ■ < a^. 

In order to compute it we follow the scheme described in Section 6 of [3]. First we define the 
function 

W(z 1 ,...,z K ;w) := (A_\T(z K ) ■ ■ ■ T(z 1 )(j) a (w)\A + ), K = d + m + m', 

where T{z) is the energy-momentum tensor as a function of the point z on the cylinder with 
the OPEs 

T(z)T(w) = ~^x'"(z -w)- 2T(w) X '(z - w) + T'(w) X (z - w) + O(l), 
T{z)<j> a {w) = -A a <j) a (w)x'(z -w) + 4>' a (w)x{z -w) + O(l) 

and 

icoth^Uy-%- 2 - 1 

where B<m are Bernulli numbers. We also need the expansion: 



1 oo 

x (z) = ±-±Y,e*', K(z)^±oo. (5.8) 



3=1 



As was discussed above we can use two different expansions for the energy-momentum tensor 
as well 
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the "local" expansion in vicinity of z = 

oo 

T{Z) = lnZ~ n ~ 2 , 
n=—oo 

the "global" expansion when $t(z) — > ±00 



T(z) = ^ - (5-9) 

71= — OO 

The action of the local Virasoro generators \ n on a local field O(w) is defined through the 
contour integral 

(inO)H= / ^(z-^r +1 r(z)OM, 

where C w encircles the point w anticlockwise. 

The conformal Ward-Takahashi identity allows to determine the function W(z\, . . . , zk) re- 
cursively: 

c K I 
W( Zl , ...,z K ;w) = -—Yl " ZjWfa, ,zk;w) (5.10) 

+ (~^x'{ z i - *j) + " z i) " " " ^X{zi - w) 



0=2 



+ (A+ - A_)x(«i - w) + ^(A+ + A_) - ^W(z 2 , ...,z K ; w). 

Of course, the term containing the difference A + — A_ drops for the case A + = A_ which is 
only interesting for us here. 

In order to calculate the above object (n^, . . . , n' m ,; a±, . . . , a^ni, . . . , n m ) given by (|5.6p . we 
proceed in several steps: 

step 1: take the recursion (|5.10p and expand x( z i ~ " • " )> x'C^i ~~ " " • ); x'^^i — ' ' ' ); x"'( z i ~ ' ' ' ) 
using the expansion (|5.8p for — >• —00 and then having in mind the expansion (|5.9p 

take there the coefficient at e ni2:i , 

step 2: repeat this procedure consequently for the variables Z2, £3, ■ ■ ■ , z m > taking every time 
the coefficients at e n 'i Z2 , e n 'z Zi , e n '™> z ™' , 

step 3: similarly we proceed with the next m variables £ m '+i, ■ ■ ■ , z m ' +m taking the expansion 
(15. 8P for 3?(z m / + i) — > 00, . . . , $t(z m ' +m ) — > 00, further using the recursion (]5.10p and picking 
up coefficients at e — ni * m '+ 1 , . . . , e _nm2m '+ m , 

step 4 : now one can easily compute the limit w — > and then apply (|5.10p with respect 
to the variable z m > +m+ i, take the local expansion (|5.7p of x( z m'+m+i)i x' { z m' +m+i) ■, 
x" ( z m'+ m +i), ■■■ and calculate the contour integral J c d^m'+m+i^"/ 2 " 1 ^^ • • • with the 
expression obtained in this way, 

step 5: repeat the step 4 with respect to the residual variables z m '+m+2i ■ ■ ■ ,Zk every time 

I Co aZm + m + 2Z m'+m+2 

^._-2a d +l . . . 

'Co 



calculating the contour integrals j Cn dz m i + m +2 z r ri'+m+2 ' ' ' e ^ c - U P ^° * ne ^ as ^ integral 

Icn d z K z K 2ad+1 ' 
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In this way we can obtain, for example, for the case d = 0, N = 2 

(2;0;2) = ^ + 4(A 2 -A + A + ), 

(1, 1; 0; 2) = (2; 0; 1, 1) = 2(A 3 - A + 3A+), 

(1, 1; 0; 1, 1) = A(A - 1)(A 2 - A + 2) + 4A+(2A 2 - 2A + 1) + 8A^, 

where again we used the shorthand notation Z\ = Aq,. 
For d = 1, ai = 1, N = 2 we get 

(2; 1; 2) = -1a(A 2 - 61A + 42) + c (-^A 2 + ^A + l) - ^ 



+ A + ( 4A 2 + ^A + ^ + 8 ) + 4A^, 



35 . c 
—A + - 
3 3 

;i,l;l;2) = (2; 1; 1, 1) = -^A(A + 1) (2A 2 + cA - 86A + 72 - 7c) 



+ A+ ( (2A + 3) ( A 2 + ^A + 4 



+ 6 A 2 



1 



(1, 1; 1; 1, 1) = A ( A - 1) ( 2A + (c - 50)A 2 - (c + 44)A + 2c - 96) 

+ — (6A 4 + 32A 3 - 2(c - 71) A 2 + 2(c + 29) A + 48 - c) 
6 ' 

A 2 

+ — t (24A 2 + 70A + 60 - c) + 8A 3 . 



(5.11) 



(5.12) 



And for d = 1, ai = 2, N = 2: 



1441 \ 



A /A 2 359 , 1921c \ 1921 AA+ 

2: 2; 2) = A — + A + 3 + -, 

V ' ' ; V 60 60 480 J 60 

(1, 1; 2; 2) = A(A + 1) (J^ + ^ A + c - l) + A (l2A 

, /A(A-l) 241 \ /121 2 599 A 1921\ A 
1, 1; 2; 1, 1 = A 2 A - 1 — i '- + + A A + A + A. 

v ' ' ' ' ' y 1 \ 240 120/ V 30 30 60 y 



40 



- 



A_ 



(5.13) 



+ 



481AA^ 
30 



6 The level M = 8 

Here we show how the procedure generally described in the previous section works for the case 
of level M = 8. As was pointed out in Introduction, there are 5 monomials of the local Virasoro 
generators which generate linear independent descendant states modulo integrals of motion. Let 
us arrange them as a vector 



Vi 



( l 4 - 2 \ 

l-4l 2 _ 2 
l 2 -4 

l- 6 l-2 

V u J 
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As for the fermionic operators (|5.4p . we also have 5 possibilities for M 
as a vector 



which we also take 



/ 



/even 
<Pl,7 

<t>if 

<t>lT 

/qCFT* fl CFT* CFT* «,CFT* 



where we used the even and odd combinations (|2.18p . We would like to determine the 5 by 5 
transformation matrix U 



y fcrmi UVh (6.1) 

where the weak equivalence "=" should be understood in the same way as in the formula (I2.2ip 
and all the matrix elements of U depend only on the conformal dimension A = A a and central 
charge c. In [5] we were not able to uniquely fix this matrix from the consideration of the ground 
state case N = 0. If we substitute the both sides of (|6.ip into (|2.ip where the matrix elements 
of U in every row correspond to the coefficients in the polynomial P a [{l-k}) and use (j5.5j) for 
k = then we get 



(A-|(CfW a (0)|A. 
(A_|0 a (O)|A + ) 



(6.2) 



where 



/ (ng + ng)/2 \ 



(nL°J + nS)/2 



nS)/(2da) 



(6.3) 



and fi^v denote the functions Qjj/(p, a) given by (|4.1ip which are calculated for the ground state 
case N = 0, k = &s was explained in Section 11 of Prom the vanishing property follows 
that all five equations in (|6.2p are polynomials with respect to p 2 of degree 4. In all 5 cases one 
of the equations is fulfilled automatically. It means that every equation leaves an one-parametric 
freedom. In other words, in every row of U one matrix element is left undetermined. 

Therefore in order to fix the matrix U completely we need to involve excitations. Here we 
restrict ourselves with the case k = 1 and N < 2 only. Let us remind the reader that the 

descendant level for the excitations N = | ^ (Ir + ' k ^ + Ir '^). So, for the case N = 1 we have 

r=l 

only one possibility = /("i 1 ) = As in the ground state case N = 0, the descendant 

space is one-dimensional because there is only one vector L_i|A+) here. When we wrote the 
paper [4] we thought that considering this excitation would help us to fix the above mentioned 
freedom. As appeared it is not the case because it does not produce any additional constraints 
on the transformation matrix U. We leave the checking of this fact as an exercise for the reader. 

So, we have to consider the case N = 2. There are two possibilities = {1}, J^' 1 ) = {3} 

and = {3}, i^"' 1 ) = {1} that correspond to the two-dimensional descendant space spanned 

by two vectors L_2|A+) and L?_ 1 |A+). Now let us consider the states (I5.ip 



|A+;{1},{3}) =AL_ 2 \A. 



+ BL 2 



-i 



A_ 



(6.4) 
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|A+; {3}, {1}) = CL_ 2 |A + > + DLi 1 |A+) I 
<A_;{1},{3}| = A(A^\L 2 + B(A„\L 2 , 
<A_; {3}, {1}| = C(A_|L 2 + 5(A_|Lf, 

where we use simpler notation for the coefficients A^ 1 ^'^ 3 ^ = A, A^^'^ 3 ^ = B, A^^'^ = C, 
^({sj.ii}) ^ Y) anc | similar for A, B, C, D. As appeared we can take them A = A, B = B, 
C = C, D = D. 

Now we should satisfy the normalization condition (j5.2[) 



<A_; {1}, {3}|0 Q (O)|A + ; {1}, {3}) = <A_; {3}, {1}|^(0)| A + ; {3}, {1}) = 1. 

Substituting the above formulae (|6.4p here, using the notation (|5.6[) and the fact that (1, 1; 0; 2) = 
(2;0;1,1), we get 

A 2 {2; 0; 2) + B 2 (l, 1; 0; 1, 1) + 2AB(1, 1; 0; 2) 

= C 2 (2; 0; 2) + L> 2 (1, 1; 0; 1, 1) + 2CD(1, 1; 0; 2) = 1, (6.5) 

where (2; 0; 2), (1, 1; 0; 1, 1), (1, 1; 0; 2) are given by (|5.1ip . Now we can use the first row 
of (|2.2ip . the formulae (|5.5p and (|6.5p in order to obtain 



(2;0;2) (1, 1; 0; 1, 1) (1,1;0;2) 

(2;1;2) (1, 1; 1; 1, 1) (1,1;1;2 

(2;2;2) (1,1; 2; 1,1) (1,1; 2; 2 

,(1,3) 




ViKi 3) -^3 3) )K 



Here we denoted ^ the function Qjj/(p,a) given by the formula (|4.1ip where 



e(ij/(2u),ij'/(2v)\p,a) 

is determined by the equation (|4.7p with k = 1 and = {1}, /C - ' 1 ) = {3} while d a is 

taken from (|2.19p . The matrix elements of the second and the third row in the left hand side 
of (|6.6|) are given by (|5.12p and (|5.13p respectively. The matrix equation (|6.6p can be solved 
with respect to A and B by inverting the 3 by 3 matrix. These equations are over deter mined. 
The first two of them give A 2 , B 2 and we can check that the last equation which gives 2^45 is 
fulfilled automatically. 

In a similar way one can get the coefficients C and D from the equation which can be obtained 
from ()6.6p by changing A — >■ C , B — >■ D and fi^jf — > ^fft where fij 3 'v is defined through the 
equations (|4"TB and (@2D with = {3}, = {1}. 

The result for few leading orders with respect to p looks as follows 

, l / ,_j c-22. ,_ 2 /A(A-l) c 2 + 148c- 860\ . ,_ 3 _ _ 4 . 

^ = -(pi/) 1 H (pi/) 2 - — H (pi/) d + Ofp 4 ), 

4 VjP 7 192 VjP ; V 16 18432 / 

5 = \(pv)~ 2 ~ ^M" 3 + 0(P~ 4 )- (6-7) 

The coefficients C and -D can be got from A and B respectively through the substitution of p 
by -p. 

As in the N = case described above, one can substitute (|6.ip into (|5.3j) and use (|5.5p in 
order to get 

A 2 {A^\L 2 (UV i )M0)L- 2 \A + ) + £ 2 (A_|L 2 ( L n/ 1 )</> a (0)L 2 _ 1 |A + ) 
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+ 2AB(A_|L?(^i)0 a (O)L_ 2 |A + ) = F (1 ' 3) , 



(6.8) 



where V^ 1 ' 3 ) can be obtained from W°) given by (|6.3|) via the substitution Qjji ^l^f 



( (^ 3) +4 1 i 3) )/2 \ 



(1,3) 



TO 



(1,3) 
3,5 



& 3) )/( M ° 

)/2 



(1,3) 
5,3 



\ (1,3) (1,3) 
V 4 1,1 Si 3,3 



/3 



)(1,3) 



Both sides of the equation (|6.8p can be represented as a power series with respect to p _1 starting 
with p°. Then one can equate coefficients standing at powers and get as many equations as 
necessary. It turned out that all the equations obtained by taking coefficients at p° up to p~ 8 
fulfill automatically. The unknown matrix elements of the transformation matrix U are fixed 
only in the order p~ 9 . All equations which stem from further orders with respect to p~ 1 should 
be fulfilled automatically. Unfortunately, so far we could not prove it or even check any further 
equations because of complexity of calculations in the intermediate stage. We hope to do it 
in future. We checked that if one substitutes A — > C, B — > D in the left hand side of the 



(1,3) 
3,j' 



fr 3 ', in the right hand side of (|6.8p then the equation obtained 



3,0 



equation f|6.8j) and fi 

is fulfilled automatically up to p~ 9 . Let us explicitly show the final result using even and odd 
combinations (|2. 18[) : 



<Pl,7 — 



1 4 „ 2 + (4(c - 4) A + 4(c + 8))/(3(A + 4))l_ 4 l 2 „ 2 

- (-43c 2 + 924c - 16340 - (21c 2 + 1102c - 18535)A 

- (11c 2 - 198c + 775)A 2 + 40(c - 25)A 3 )/(45(A + 4)(A + ll))l 2 



12(c 2 



- 1540c + 17264) - 2(43c 2 - 298c + 2652)A 

31c + 574)A 2 + 32(c - 28)A 3 )/(l5(A + 4)(A + ll))l_ 6 l_ 



(-45c 3 + 1637c 2 - 137176c + 2033360 - 20(2c 3 + 261c 2 - 7623c + 67411)A 
4(10c 3 - 439c 2 + 7142c - 5825)A 2 



/odd rsj 

<Pl,7 — 



+ 96(2c 2 - 81c + 705)A 3 )/(105(A + 4)(A + 11))1_ 8 , 
(4A)/(A + 4)1_ 4 1 2 _ 2 
+ (20(c - 28) + 24(c + 47)A + 8(2c - 11)A 2 )/(15(A + 4)(A + ll))l^ 4 
+ (-40(c - 28) + 8(13c - 44)A + 16(c - 8)A 2 )/(5(A + 4)(A + ll))l_ 6 l_ 2 

- (20(c - 28) (c + 79) - 4(15c 2 + 1316c - 20143)A 

- 4(15c 2 - 398c + 4359)A 2 + 32(c - 33)A 3 )/(35(A + 4)(A + 11))1_ 8 , 
= 1-2 + (8(2c + 13) + 12(c - 16)A)/(9(A + 4))l_ 4 l 2 _ 2 

+ (2(262c 2 - 4271c + 82750) + 3(59c 2 + 2338c - 67785)A 

+ (79c 2 - 1502c - 6965)A 2 + 120(c - 25)A 3 )/(405(A + 4)(A + ll))l 2 4 

+ (4(68c 2 + 7571c - 86380) + 4(188c 2 - 3379c + 10310) A 

+ 24(4c 2 - 127c + 1125)A 2 )/(135(A + 4)(A + 11))1_ 6 1_ 2 

+ (2(420c 3 - 9013c 2 + 711929c - 10449400) 

+ 4(70c 3 + 9741c 2 - 370938c + 3325745)A 



(6.9) 
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+ 20(14c 3 - 629c 2 + 9532c - 48805)A 2 

+ 96(c - 25)(2c - 31)A 3 )/(945(A + 4) (A + 11))1_ 8 , (6.10) 
< 5 d -(4A)/(3(A + 4))l_ 4 l 2 _ 2 

+ (-4(17c - 236) - 24(c - 43)A + 8(4c - 127)A 2 )/(l35(A + 4)(A + ll))^ 
+ (8(17c - 236) + 8(47c - 536)A + 48(c - 18)A 2 )/(45(A + 4)(A + ll))l_ 6 l_ a 
+ (28(15c 2 - 571c + 7708) + 4(35c 2 + 2708c - 66859)A 
+ 20(7c 2 - 222c + 767)A 2 + 96(c - 33)A 3 )/(315(A + 4)(A + ll))l_ 8 

and for the fourth-order combination we get 

^CFT^CFT^CFT^CFT* ^ 1/12 . ^ + _ 54 + 2(c _ 22 )A)/(l8(A + 4))L_ 4 1 2 _ 2 

- (43c 2 - 1426c + 11664 + 2(15c 2 - 302c - 1069)A 

+ 2(c 2 + 86c - 2667)A 2 + 16(c - 25)A 3 )/(216(A + 4)(A + 11)) £ 4 

+ (215c 2 - 8714c + 74952 + (125c 2 - 3856c + 13208)A 

+ 2(5c 2 - 73c - 1816)A 2 + 16(c - 28)A 3 )/(l80(A + 4)(A + ll))l_ 6 l_ 2 

- (-25c 3 + 1755c 2 - 39410c + 326592 - 4(30c 2 - 1393c + 9520)A 

+ 8(llc + 289)A 2 + 192A 3 )/(72(A + 4)(A + ll))l_ 8 . (6.11) 

The elements of the above matrix U can be easily got from these data. It is interesting to note 
that the determinant of U has relatively simple factorized form 

det U = -32(c - 25)(c - 28) (c - 33)(A - 1) (6.12) 
(c + 2 - 2(c + 11)A + 48A 2 )(-c 2 + 18c + 175 + 16(c - 25)A + 192A 2 ) 
X 382725(A + 4)(A + 11) ' 

The numerator of (16. 12ft is proportional to the following product 

{y - 2) 2 {y - 3)(i/ - 4)(2i/ - 3)(3i/ - 4)(ai/ - l)(av - 2)(au + 1 - v)(av -\-v) 

x(av + 2- v)(av - 2 + v)(av - 2 - v)(av + 1 - 2v)(av + 2- 2v)(av + 2- 3u). 

It would be interesting to understand the meaning of the degeneration points like ^ = 2,^ = 3, 
v = 4, 1/ = 2/3, z^ = 3/4 and a = 1/u, a = 2/v, a = — (1 — v)/v etc. where the determinant 
det U = 0. 



7 Conclusion 

In this paper we demonstrated that the method developed in [1] works for excited states as well. 
As we saw with the example of the descendant level M = 8, it is not possible to completely 
determine the transformation matrix between the usual basis constructed through the action of 
Virasoro generators on the primary field and the fermionic basis constructed in [lj [2] without 
involving excitations. On the other hand the matrix elements of the transformation matrix 
should not depend on the fact which excitation is taken. This should provide an interesting 
compatibility condition on the structure of the three-point functions both from the point of 
view of the CFT and the lattice XXZ model we started with. In this paper we were able to 
treat only the case of excitations corresponding to the descendant level N < 2. It would be 
interesting to check the above mentioned compatibility for the case of other excitations with 
N > 2 and for the higher descendants M > 8 as well. Of course, it is very important to find 
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a general proof of the compatibility condition. We think it is also interesting to study singular 
points of the transformation matrix mentioned in the end of the previous section. This may 
shed new light on the structure of the Virasoro algebra, Verma modules and singular vectors. 
One more important generalization of the results obtained here, which is still out of our reach, 
would be to treat the case of general values a, k, k' and also the case of different excitations 
inserted at +00 and —00 on the cylinder. In both cases the function p is not 1 and we would 
have to generalize the whole Wiener-Hopf factorization technique. We hope to return to these 
questions in future publications. 

A The function <3>( p ) and integrals of motion 

Here we show several further orders of the 1/p-expansion of the functions used in Section [3] 
for the case of excitations with k = 1. For simplicity we use shorter notation mo = 
m\ = l[ . First let us show expressions for several coefficients ^ (Z) of the expansion fl3.26f) 
for the function 

$Wm = - 

W ° [l) l(l + i/(2u)Y 

*! p) (o = o, 

^ J \l) = ^(-l + l2(m + m 1 )), 
*P(l) = - l -^^(rnl-rnl), 

tf^CO = - 2 ^(i (2 -l) - v){V2 + 60(m ° + m?)) + 5(21/2 - 11(1 - u))im ° + m?) 

+ 5(1 + u)(2u - l)(mo + mi)(6(mo + m\) - 1) + v 2 + 3(1 - v)), 

^ = 147w|l-,) 2 (m ° " -?)( 384 ^- 2 ( 1 " -) 2 K + mj) 

- 16ilv(l - v)(b(2v 2 + 11(1 - v))(ml + m 2 ) 
+ 2(1 + v)(2v - l)(12(mo + mi) - 1)) 

- (20z/ - 220z^ 3 + 68I1/ 2 - 922^ + 461) (m 2 , + m?) 

- 2(1 + v)(2v - \){2u 2 + 23(1 - z/))(12(m + mi) - 1)) 

*i p) m = <^640i/ 3 (-31 + 252(77^ + m?))(l-z;)V 

6 y> 8360755200^ 3 (1 - vf 1 v y ljn ' 

+ 1440Z 2 (1 - v) 2 v 2 (b§7{2v 2 + 11(1 - v)){ml + m\) 

+ 21(1 + u)(2u - l)(m + mi) (7 - 10 (m§ + ?n 2 - m mi) + 120(mo + m?)) 

- 164z^ 2 - 755(1 - u)) - 90Uv(l - v)( 21(6085 - 12170i/ + 8769i/ 2 - 2684z^ 3 



+ 244i/ 4 )(m§ + ml) + (1 + i/)(2i/ - l)(420(34z/ 2 + 295(1 - i/))(mg + mf) 
+ 840(14i/ 2 + 113(1 - v))(mlmi + mom 3 ) 

- 2520(2i/ 2 + 23(1 - v))mlml + 70(274i/ 2 - 257(1 - v))(m? Q + mf ) 

+ 2520(1 + v)(2v - l)(m + mi)(-m - mi + 12m mi) 

+ 7(158z^ 2 + 761(1 - 1/)) (mo + mi)) - 96(159 - 318i/ + 240i/ 2 - 81z^ 3 + 8v 4 ) 
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- 21(210403 - 6312091/ + 813759z> 2 - 575503z/ 3 + 219354z/ - 36804i/ 5 
+ 1288u e )(m 5 + ml) - (1 + u)(2u - 1) (420(20281 - 40562^ + 25053z> 2 

- 4772i/ 3 + 52z/ 4 )(m 4 + m\) - 840(-6407 + 12814i/ - 6771i/ 2 + 364i/ 3 
+ 436z/ 4 )m mi(mg + m\) - 2520(2489 - 4978i/ + 3837i/ 2 - 1348z/ 3 

+ 308z/ 4 ) mgm? + 70(-31607 + 63214i/ + 16989z/ 2 - 48596i/ 3 

+ 3316i/ 4 )(?ng + m\) + 12600(1 + v)(-l + 2v){2u 2 + 35(1 - v)){m + mi) 

x (-mo - mi + 12m mi) - 7(-39599 + 79198i/ - 52347z/ 2 + 12748i/ 3 

+ 2452i/ 4 )(?n + mi)) + 103680(1 - i/)(5 - 10^ + 10^ 2 - 5z/ 3 + i/ 4 ) I. (A.l) 



One can check that the ^-function in the ground state case is reproduced if one takes mo = 
mi = 0. There is a connection of the function with the integrals of motion described in [7]. 
In the ground state case this connection was given by the formula (10.17) of [3J: 



'2n-l 



^(P) f i( - 2n ^ *M n (2n - l)(2 l / 2 )»-y n - 1 . (A.2) 



In the case of excitations the integrals of motion are given by matrices. Let us consider the first 
three integrals of motion I±, 1$, 1$. Their explicit expressions via the Virasoro generators can 
be found in [6]. 

Again let us consider an example of the excitations with N = 2 and k = 1 with two pos- 
sibilities = {1}, ^ = {3} and = {3}, = {1} which correspond to the 
two-dimensional descendant space spanned by two vectors L_2|A + ) and L 2 1 |A + ) or their linear 
combinations | A + ; {1}, {3}), |A + ;{3},{1}) given by the formula (|6.4D . 

Let us start with the simplest case of the very first integral of motion I\ given by the 
formula (11) of the paper [Bj: 

U 24 

This operator is diagonal for the above two-dimensional space 



L_ 2 |A+)\ _ / cx /L_ 2 |A, 

L 2 _ 1 |A + )J-l A+ + 2 -24J [l^A, 



and 



c 47 

A+ + 2 = 2{vvf A . 

+ 24 Ky ' 24 

We can easily check that we get exactly the same result if we substitute the expansion (|3.26p 
with the coefficients given by (|A.lj) for / = i/{2v) and mo = 1, mi = 3 or mo = 3, mi = 1 into 
the formula (|A.2|) with n = 1. We see that in this case all ^\i/(2u)) = with j > 3 even 
without fixing mo and mi. 

The next case of 1$ is a bit less trivial since ^3 is not diagonal any longer. Following [6J, we 
have 



1% — 2 \ L— n Ln A~ Ln — Ln A — 

^ 12 2880 



n=l 
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We see that only the first term here is not diagonal. Using the Virasoro algebra, one can check 
that 

L_ 2 |A+)\ f c 6\ /L_ 2 |A H 
L 2 _ 1 \A+)J ^12A + 4j {l^A, 

+ f 8A+ + (A + + 2? - i^±M^±±3 + 4* + 22 ) A f^ 2 \A 



12 2880 J \L 2 _i\A. 



We can diagonalize the matrix 



c 6\ (C D\ 1 f\ + \ (C D 



12A+ A) \A B) V0 A - / V^4 -B 
with the eigenvalues 



6t/ 2 _ 5(1 _ ^ 3^/64^(1 - gg + (1 + ^) 2 (2^ - l) 2 
A " 2(1-!/) 2(1^) (A '^ 

and the matrix elements A, B, C, D taken from the definition of the vectors | A + ; {1}, {3}), 
|A + ;{3},{1}) given by (|6.4p . Hence these vectors are eigenvectors of the matrix I3. Several 
leading terms of 1/p-expansion of A, B where shown in (16.71) w hile C, D can be got from A, B 
by changing p — > —p. So, for the eigenvalues of ^3 we obtain 



J 10 



/« = A± + 8A + + (A + + 2? - iE±m±±V + ^ 
= 4M* + ^M^12 P , «04^- 5769(1 -,) 



2 Kl ' * 960(1 - v) 

3 (l + z^) 2 (2zy- l) 2 . , 1 _ _ 2 , .. 
^32 (1-zv) 2 {W) )■ ( } 

We can check at least for number of first orders that we get the same 1/p-expansion if, like 
in the previous case we substitute the expansion (|3.26|) this time for I = 2>i/{2v) and tuq = 3, 
mi = 1 in case of Jg and mo = 1, mi = 3 in case of 1% ^ into the formula (1A.2|) with n = 2. 
In contrast to the case of the first integral of motion I\, the series expansion (|A.4p does not 
terminate because of the square root in the expression (1A.3|) for the eigenvalues X^. Also we get 
both even and odd powers with respect to p in contrast to the case N = corresponding to the 
ground state. 

The next case can be treated similarly. First we take 1$ from [6] 



/c + 11 2 c \ 
'■ L ni L n2 L n3 : + I — g — n — 1 — ^ ) L^ n L n 

ni+ri2+n%=0 n=l ^ ' 



3^ r c + 4 r2 (c + 2)(3c + 20) c(3c + 14) (7c + 68) 

2 ^ 8 576 290304 

r=l 

with the normal ordering : : for which the Virasoro generators with bigger indices are placed to 
the right. Then using the Virasoro algebra, one can come to the following formula 

i 2 -i|A + )J - 5 I 6 + 24 + A+ ) 1 12A+ A) U^IA. 



10 We hope the reader will not mix these eigenvalues /■ with the integrals of motion mentioned in Section[2] 
that equal to each other when k' = k and A+ = A_ . 
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236 -c a2 / c 2 227c 3845 \ . 

+ Ai H Ai + 1 A i 

V 8 V 192 288 72 y H 

/ 3 c 3 1361c 2 _ 7325c 221 \ /L_ 2 |A„ 

+ V + ~ 13824 + 145152 ~ 5184 + ~36~ ) V L -il A - 

We see that the matrix I5 can be diagonalized by the same similarity transformation as I3 
because all the integrals of motion commute. The result for the asymptotic expansion of the 
two corresponding eigenvalues I 5 looks 

r(± ) , , 6 235. , 4 1(W , 3 4808i^ 2 - 12503(1 - v) . 

H } - 8(H 6 + — [pvf T 120(p^) 3 ^ _ -A > (p„)2 

15(4i/ 4 + 36z/ 3 + 21^ 2 - lUv + 57) 
16(1 - uf 



=F T7^ 7T\2 (P") 



822432i/ 4 + 2338220^ 3 - 821605^ 2 - 3033230z^ + 1516615 . _ n 
+ 96768(1 - ,) 2 + °^")- 

Again this result perfectly matches the expansion of (lA~2l) for n = 3 with tfW (5i/(2i/),p) taken 
for mo = 3, mi = 1 in case of Ig and for mo = 1, mi = 3 in case of Jg . Similar to I3 , the 
integral I5 and all further integrals I^n-i are no * polynomials with respect to j> in contrast to 
the ground state case. 



B The functions F(x,p) and x^(p) 

In Section [3] we discussed several functions defined within the TBA approach. In the previous 
appendix we showed few orders of the function ty(p\l,p). We also need the functions F(x,p) 
and x^r(p). As was described in Section El the asymptotic expansions of these functions and the 
function ^^ p '(l,p) are calculated order by order with respect to 1/p via the iterative procedure. 
There we explained how does the very first iteration work. Once the function ^f^ p \l,p) is found 
up to some order, the next order of the function F(x,p) can be obtained via the equation f|3. 1T|) . 
The coefficients of (|3.18j) are determined from the equation fj3.21f) . The further iteration steps are 
straightforward but the answer becomes rather cumbersome already after several iterations. For 
the reader who wants to check his own calculations we show few orders of the expansions (13.231) 
and (|3.18p in the case k = 1. We do not think it would be instructive to show further orders. 
As in the previous appendix we use the shorthand notation mo = Ii , mi = l[ • 

- /x (1 + i/)(2i/-1), , . . 9 2i/ 2 + 11(1 - u) 

^ = 57M1-.) " - 12(m » + - " 48,(1-,) ■ 

_ ,(1 + ^-1)^ + 23(1-,)) _ ^(1 + ^(2,-1); _ + 
u ' 9216z/ (1 — is) 27648zy 2 (l - v) 2 K K u n 

3 4z/ - 44^ 3 + 309z> 2 - 530z^ + 265 
+ X 6912zy 2 (l-^) 2 ' 

h{x) = - 23^78^V-W ( 12 ° (436,/4 + 364 ^ " 67711/2 + U8U " ~ 64 ° 7) 
x (mjj + m\) + 3600(1 + v)(2u - l)(2v 2 + 23(1 - u))(m + mi) 

" (1 - 6(m + mi)) + 2452i/ 4 + 9U8u 3 - 50547^ 2 + 82798i/ - 41399^ 

a(l + u)(2u - 1) 4 _ 1676zy3 + 2859^ 2 - 2366z^ + 1183) (ml - m? 

3317760i/ 3 (l - ^) 3V A 



2 2\ 

, -m x ) 
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19906560^ 3 (1 - v) 



(556z/ 4 - 2036^ + 3039z/ 2 - 2006^ + 1003) 



• 4 

IX 



ix 2 (l + u)(2v-l) fKKRtA on ^_ 3 

x (1 - 12(m + mi)) ^7 -^(1112i/ 6 - 2796i/ 5 - 5154z/ 4 

v v. u-r ljj 99532801/3(1 _ j/)3 v 

+ 64603z/ 3 - 154059z/ 2 + 146109i/ - 48703), 
(l + v )(2v-l) ,_ 2 
15288238080i/ 4 (l - i/) 
+ 233059i/ 3 - 658107i/ 2 + 649893z/ - 216631) (mg + m\) 
- 4(1 + u){2u - l)(932i/ 4 + 788z/ 3 - 22227^ 2 + 42878i/ - 21439) 



F 3 (x) = -i i. 1 + v )( 2v jj / m 2 _ 2) C 3 (773 6l/ 6 _ 49644^5 + 41430l/ 4 

dV 7 15288238080i/ 4 1 - v 4 V 17 V v 



x (1 - 12(m + mi)) ) + x ( l + u )( 2u jj — / 10 80(1496i/ 6 - 12924i/ 5 

V V u-r lyyy -r 2 29323571200i/ 4 (1 - z/) 4 V v 

+ 42270z/ 4 - 75161i/ 3 + 78753^ 2 - 49407z/ + 16469) (m 3 , + m\) 

- 960(1 + v)(2v - l)(556z/ - 1676z> 3 + 2859i/ 2 - 2366z/ + 1183) (mo + mi] 

x (1 - 6(m + mi)) + 138728i/ 6 - 926052i/ 5 + 2802450z/ 4 - 4743023z/ 3 

9 \ 9 fl + z/)(2z/ - 1) , 9 9x 
+ 4847079// 2 - 2970681^ + 990227 + ix 2 — - ; -- r - '— (ml - m?) 



637009920i/ 4 (l - z/) 4 
x (4568z/ 6 - 40572z/ 5 + 138270z/ 4 - 239513z/ 3 + 230049z/ 2 - 132351z^ + 44117) 



l + v) 2 (2u l) 2 , or>Q/|Tj4 cn sM,,3 , 10111, .2 
5733089280i/ 4 (l - z/) 4 



+ 1 rwononoon.4/1 . .u ( 2284t/ ~ 8084i/ 3 + I2111i/ 2 - 8054z^ + 4027) 

n5 



x (1 - 12(m + mi)) + -(9136z/ 8 - 76576i/'' + 196840z/ B 

v v i;y 4777574400z^ 4 (l - i/) 4 ^ 

+ 149096i/ 5 - 3382325z/ 4 + 10503740z/ 3 - 14391926i/ 2 + 9334868z/ - 2333717), 
F A (x) = j} 1 + ^} 2U ~ yj 7F f 6048(4048z/ + 208064i/ 7 - 1145688i/ 6 + 1757040zv 5 



2157476157849600z/ 5 (l - vf . 
+ 1432457z/ 4 - 8416452z^ 3 + 12210950z/ 2 - 8061828z> + 2015457) (m^ + m\) 
+ 3809872z/ 8 - 215990816i/ 7 + 939983096z/ 6 - 1142084632z/ 5 - 1598906179z^ 4 
+ 6538494484i/ 3 - 8737091978i/ 2 + 5620794700^ - 1405198675 

- 196(1 + u)(2u - l)(m +mi)(240(11528i/ 6 - 109308i/ 5 + 124842z/ 4 + 705307z^ 3 

- 2193591z^ 2 + 22091251/ - 736375) (m 3 , + m\) + 120(27112i/ 6 - 287566z/ 5 

+ 627666z/ + 187193z/ 3 - 2262129z/ 2 + 2602239^ - 867413)m mi(m 2 1 + m\) 
+ 10(247064z/ 6 + 795588i/ 5 - 9542418z/ 4 + 12330151i/ 3 + 6743697i/ 2 - 15490527^ 
+ 5163509) (m 2 , + m\) + 10(575464i/ 6 + 728508z/ 5 - 17201838z/ 4 + 25221881i/ 3 
+ 6701007i/ 2 - 23174337i/ + 7724779)m mi + 360(2i/ 2 + 23(1 - i/))(676i/ 4 

- 3716z/ 3 + 12489z; 2 - 17546i/ + 8773)m mi(m + mi) - 5040(1 + v){2v - 1) 
x (68z/ 4 + 92z^ 3 - 2223i/ 2 + 4262i/ - 2131) (m + mi) + 7(35272i/ 6 

- 272436z/ 5 + 362346z/ 4 + 559733z/ 3 - 2128749i/ 2 + 2218659^ - 739553)) 



x(l + v)(2v — 1) . 9 0s / . fi 7 R 

— i ^ — — ^{ml - m? 8 114704z/ 8 - 186848i/ 7 - 1703896z/ 6 

1284211998720i/ 5 (l - z/) 5V 17 V v 

+ 8391976z/ 5 - 17873683z/ 4 + 22437364i/ 3 - 17519338z/ 2 + 8605900i/ - 2151475) 
x (ml + m\) + 7(1 + u){2u - l)(11240z/ 6 - 137412i/ 5 + 480594i/ 4 - 893735i/ 3 
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965295z/ 2 - 6221131/ + 207371)(1 - 12(m + mi)) 
x 2 (l + u) 2 (2u- l) 2 



i 



(120(1247080z/ 6 - 6770580z/ 5 + 15390258z/ 4 



77052719923200z/ 5 (l - vf \ 

- 18566683z/ 3 + 12601659z/ 2 - 3981981z/ + 1327327) (m 3 , + m?) 

+ 1680(4568z/ 6 - 40572z/ 5 + 138270z/ 4 - 239513z/ 3 + 230049z/ 2 - 132351z/ 

+ 44117)(m + mi)(l - 6(m + mi)) + 7(1201400z/ 6 - 6364860z/ 5 + 14007558z/ 4 

- 16171553z/ 3 + 10301169z/ 2 - 2658471z/ + 886157)) 

x 3 (\ + v\{1v — \\ 

+ \ -r^(ml - m 2 )(2622064z/ 8 - 14133568z/ 7 + 33585976z/ 6 

321052999680z/ 5 (l - z/) sV 1A 

- 49220560z/ 5 + 55866331z/ 4 - 52813996z/ 3 + 43115890z/ 2 - 21866764z/ + 5466691) 

ix A (l + u)(2u - 1) o 7 fi 

i ^Hr7 — ^— ^ 2622064z/ 8 - 14517280z/ 7 + 35135944z/ 6 

3852635996160z/ 5 (l -z/)5 l 

- 51550552z/ 5 + 56686759z/ 4 - 51453364z/ 3 + 41083762z/ 2 - 20513692z/ + 5128423) 



x (1 - 12(m + mi)) — tt(5244128i/ 1 - 25645072z/ 

v v u-r uj 4815794995200i/ 5 (l - z/) 5 ^ 

+ 45932736i/ 8 - 30480936z/ 7 + 22676526z/ 6 - 291308781i/ 5 + 114045686: 

- 2039726202z/ 3 + 1941710088z/ 2 - 954519025z/ + 190903805). 



6 



9 



oo 

For the parameters xf(p) = x^(p) = ^ x^p~i which determine the Bethe roots via the rela- 
tion (|3.18p we had the initial conditions xj" = ^rp, x^ = — ^f 1 and then 

= ,(2^ + 11(1 ,)) _ id + ^-l) _ 

1 192i/(1-i/) 576^(1 - 1/) v v 77 

+ _ i(20z/ 4 - 220z/ 3 + 681z/ 2 - 922z/ + 461)mg 

X<2 ~ 55296z/ 2 (l - z/) 2 

i(l + i/)(2i/-l)(2i/ 2 + 23(l-i/)) / 9 o x 

+ ,, 9Qfi a n ^ -mo + 18m 2 + 12m mi - 6m? , 

55296z/ z (l — vy 



i(1288i/ 6 - 36804z/5 + 2193541/ 4 - 575503z/ 3 + 813759z/ 2 - 631209z/ + 210403)m 

159252480i/ 3 (l - z/) 3 
i(l + i/)(2i/ - l)(52z/ 4 - 4772z/ 3 + 25053z/ 2 - 40562z/ + 20281)m 3 



+ 



+ 



9953280z/ 3 (l - z/) 3 

i(l + i/)(2i/ - l)(436z/ 4 + 364i/ 3 - 6771z/ 2 + 12814z/ - 6407)mi(3mg + m 2 ) 

19906560z/ 3 (l - z/) 3 
i(l + z/) 2 (2z/ - l) 2 (2z/ 2 + 35(1 - i/))(-mo - mi + 12m mi + 6m?) 



663552i/ 3 (l - i/) 3 

i(\ + v\(2v- lV3316z/ 4 



i(l + z/)(2z/ - l)(3316z/ 4 - 48596z/ 3 + 16989z/ 2 + 63214z/ - 31607)?n 2 

79626240z/ 3 (l - vf 
i{\ + v)(2v - l)(308z/ 4 - 1348z/ 3 + 3837z/ 2 - 4978z/ + 2489)m m 2 

3317760z/ 3 (l - z/) 3 
i(l + z/)(2z/ - l)(2452z/ 4 + 12748z/ 3 - 52347z/ 2 + 79198z/ - 39599) 
2388787200z/ 3 (l - uf ' 
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>' ~ 458647142400,4 (1 -,)4 V- 3 ( 18704 - 8 + 1330336 ^ " 16442440 ^ + 75087544 ^ 

- 183046975z/ + 267969460i/ 3 - 244932514i/ 2 + 133379452i/ - 33344863) m| 

- 150(1 + v)(2v - l)(3256i/ 6 + 51156i/ 5 - 738570z/ 4 + 2884979i/ 3 - 5217867i/ 2 
+ 4530453z/ - 1510151) m 4 + 30(1 + v)(2v - l)(53816i/ 6 - 1086444i/ 5 + 6642390^ 

- 7848221z/ 3 - 4235067z/ 2 + 9791013i/ - 3263671) ml - 90(1 + v)(2v - l)(7736z/ 6 

- 49644z/ 5 + 41430i/ 4 + 233059z/ 3 - 658107z/ 2 + 649893z/ - 216631) mi (4m 3 , - m\] 
+ 7200(1 + vf(2v - lf(2v 2 + 23(1 - v)) 2 (-3m 2 - 2m m 1 + 36mQmi + m\ 

- 12m?) + (1 + i/)(2z/ - l)(74632z/ 6 - 808308z/ 5 + 3707610i/ 4 - 8492107z/ 3 
+ 10979811i/ 2 - 80805091/ + 2693503)m - 60(1 + v){2u - l)(9976z/ 6 - 100044z/ 5 

- 431430i/ 4 - 1092901z/ 3 + 1621773z/ 2 - 1290387i/ + 430129) m m 2 (3m - 2m x ) 



+ 

The coefficients xj can be got from x^ obtained from x+ by the replacement mo mi as 
follows: 

x 2j — x 2j' 2j+l ~ x 2j+l' J u, 1, . . . . 

C The function u; sc 

In Section |4] we explained how we determine our main object: the function u sc . The asymptotic 
expansion for co sc is related to the function G by the formula (|4.9p . In it's turn the coefficients Q n 
of asymptotic expansion (14. 6p of the function 



m\p, a) ~ Qj(^) m\a)p J 
i=o 

can be found with help of the TBA data of Section [3] and Appendices [Aj |B] via solving the 
equation fj4.7|) by iterations. In this appendix we show the result for several leading coefficients 
again in the case k = 1. With the exception of the coefficient Go, all other coefficients Qj are 
polynomials with respect to I and m: 

i 

®o(l,m\a) = -— , 

I + m 

Qi(Z, m\a) = 0, 

^„ . N v 2 a(2 - a) + 2i(l - u)(-i + 2v(l + m)) , 

e 2 (Z,m q) = i '- n ; n J \ i 11 1 - 12 (m + mi)), 

9o^(l — v) 

2 2 

e 3 (Z, m|a) = - ci r° 9 ,T m \ 9 f-3^V(2 - a) 2 - 24*a(2 - a)(l - z/)z/ 2 (-i + + m )) 
51zz/ z (l — V 

+ 32(1 - vf(-i + u(l + m)) (-t + 2i/(Z + m))) , 

4 (Z, m|a) = - 184320t/ 3( 1 _ 7X3 (~5^ 6 (5(m 3 + m 3 ) - 42(m + mi) 2 + 7(m + mi)) 

x a 3 (2 - a) 3 + 10iz/(l - i/)(-24(-3i + 2(1 + m)i/)(mg + m\) 

+ 2?n (0 ' 1) (-16i + 9(/ + m)u) - (-i + (I + m»)a 2 (2 - a) 2 + 8u(l - v) 

x (2u(l - u)(7 + 120(mjj + m?))(-« + (Z + m)i/) 2 - 3i/(l - i/)(l + 10m (0,1) ) 

x (-i + 2lv)(-i + 2mi/) - 20ii/(l + v)(2v - l)(-i + (I + m)i/)m (0 ' 1) 
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- iv(-i + + m)u){Av 2 + 5(1 + 24(mjj + m?))(l - v)))a(2 - a) 

- 16^ 4 (1 - z/)(2 - u)(l + 10m (0,1 ))(J - m)a(l - a)(2 - a) + 8i(l - ^) 3 
x (7 + 120(mjj + ro?))(-i + (I + m)v)(-i + 2(/ + m)v)(-3i + 2(/ + m)v) 

+ 8(1 - z^) 2 (l + v)(2v - 1)(1 + KW ' 1 )) (2(-t + (Z + m)v)(-i + 2(1 + m)u) 
+ (-i + 2lv){-i + 2mu))\ 

where m/ ' 1 ) = + m\ + 6(mo + mi) 2 — mo — mi. Unfortunately, we cannot show all the 
orders in the expansion of and other functions from the previous Appendices that we needed 
in order to get the formulae (|6,9p - (|6.1ip because the expressions become very cumbersome with 
growing order. Again the ground state result can be obtained by taking tjiq = m\ = 0. 
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